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Abstract:  Combining  the  Mie  scattering  theory  and  a  transfer  matrix 
method,  we  investigate  in  detail  the  scattering  of  light  by  spherical  Bragg 
“onion”  resonators.  We  classify  the  resonator  modes  into  two  classes,  the 
core  modes  that  are  confined  by  Bragg  reflection,  and  the  cladding  modes 
that  are  confined  by  total  internal  reflection.  We  demonstrate  that  these  two 
types  of  modes  lead  to  significantly  different  scattering  behaviors. 
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1.  Introduction 

Dielectric  microspheres,  which  can  support  whispering-gallery  modes  with  quality  factors  as 
high  as  109,  have  been  extensively  investigated  since  the  1970’s  [1-5].  Recently,  they  have 
been  proposed  for  many  potentially  important  applications  in  wavelength  division 
multiplexing  (WDM)  (as  add-drop  filters)  [6,7],  and  in  cavity  quantum  electrodynamics 
(QED)  [8,9].  A  common  tool  in  characterizing  properties  of  the  dielectric  microspheres  is 
through  the  measurement  of  light  scattering  by  the  microspheres.  Typically,  the  scattering 
spectrum  exhibits  many  sharp,  well-defined  peaks,  which  can  be  traced  to  the  presence  of  the 
high  Q  whispering  gallery  modes.  From  the  position  and  the  width  of  such  spectral  peaks,  we 
can  extract  the  resonant  frequencies  and  quality  factors  of  the  corresponding  whispering- 
gallery  modes. 

Besides  total  internal  reflection,  it  is  also  possible  to  construct  a  multi-layer  dielectric 
sphere  and  use  Bragg  reflection  to  confine  light  in  the  spherical  cavity.  In  the  literature,  there 
have  been  a  few  theoretical  discussions  about  spherically  symmetric  Bragg  resonators  [10- 
12].  In  a  recent  publication,  we  proposed  to  approximate  the  completely  spherically 
symmetric  Bragg  resonator  with  an  onion-shaped  spherical  Bragg  resonator  [13,14].  The 
onion  resonator  has  been  realized  at  Sandia  National  Laboratories  through  a  fabrication 
procedure  combining  etching  and  chemical  vapor  deposition  (CVD).  A  scanning  electron 
microscope  (SEM)  image  of  the  onion  resonator  is  shown  in  Fig.  1. 
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Fig.  1 .  A  SEM  image  of  the  onion  resonator. 


As  in  the  case  of  dielectric  microspheres,  through  investigating  the  scattering  of  the  Bragg 
onion  resonators,  we  can  infer  useful  information  about  the  properties  of  the  onion  resonator 
modes.  In  Ref.  [14],  we  have  shown  that  onion  resonator  modes  can  in  general  be  classified 
into  two  types,  the  Fabry-Perot  type  core  modes,  and  the  whispering-gallery  type  cladding 
modes.  The  core  modes  are  confined  by  the  cladding  Bragg  reflection  and  are  mostly 
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concentrated  within  the  onion  core  region.  The  confinement  mechanism  for  the  cladding 
modes,  on  the  other  hand,  is  total  internal  reflection.  Consequently,  the  cladding  modes  are 
mostly  distributed  within  the  cladding  layers  of  the  onion  resonators.  Later  in  this  paper,  we 
will  demonstrate  that  the  scattering  resonances  of  these  two  different  types  of  modes  have 
distinctively  different  scattering  characteristics. 

The  paper  is  organized  as  follows:  In  section  II,  we  present  a  transfer  matrix  theory  for  the 
calculation  of  the  scattering  spectrum  of  Bragg  onion  resonators  (assuming  a  plane  wave  as 
the  incident  wave).  With  a  slight  modification,  the  transfer  matrix  formalism  can  also  be 
applied  to  the  more  realistic  cases  where  the  incident  wave  has  a  Gaussian  beam  profile  (as  in 
most  optical  scattering  measurements).  In  section  III,  we  apply  the  transfer  matrix  theory  and 
calculate  the  scattering  spectra  of  the  Bragg  onion  resonators.  We  discuss  various  significant 
features  of  the  Bragg  onion  resonator  scattering  spectra.  We  conclude  this  paper  in  section  IV. 

2.  Theory  of  light  scattering  by  a  multilayer  sphere 


air 


Fig.  2.  Light  scattering  by  a  spherical  Bragg  “onion”  resonator. 

In  order  to  obtain  analytical  formulas  for  the  scattering  parameters  of  the  onion  resonator,  we 
ignore  the  presence  of  onion  stem  (as  shown  in  Fig.  1)  and  assume  that  the  onion  stmcture 
possesses  complete  spherical  symmetry.  Under  this  assumption,  we  can  label  the  onion 
resonator  mode  according  to  the  angular  quantum  number  l  [13,14],  Furthermore,  we  can 
classify  the  onion  resonator  modes  as  either  TE  polarized  or  TM  polarized,  where  the 
electrical  field  of  the  TE  mode  and  the  magnetic  field  of  the  TM  mode  only  have  transverse 
components  (along  the  ee  and  e,p  direction). 

The  Mie  scattering  theory,  which  was  developed  over  almost  a  hundred  years  ago  [15], 
has  been  used  extensively  for  the  scattering  analysis  of  dielectric  microspheres.  After 
incorporating  some  results  from  a  transfer  matrix  theory,  we  can  apply  directly  the  original 
Mie  scattering  theory  to  the  case  of  Bragg  onion  resonators.  We  begin  our  analysis  by 
assuming  that  the  incident  wave  is  a  plane  wave  at  frequency  co/2 n  with  unit  amplitude  and 
polarized  along  the  x  direction.  The  incident  wave  can  be  expanded  in  terms  of  various 
multipole  moments: [16] 
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where  Xl  m  =  l/i  ■  Fx  VT^h(<9,  <j>)  /  s]l(l  +  l)  is  the  normalized  vectorial  spherical 
harmonic  function  (the  angular  parameters  0  and  (|)  are  defined  in  Fig.  2),  k  =  col  c  is  the  wave 
vector  and  jj  ( kr )  is  the  /th  order  spherical  Bessel  function.  Without  a  loss  of  generality,  the 
scattered  field  can  also  be  expanded  in  the  form  of: 


ESc  =  X' 

1=1 


./  4^(21  + 1) 


^hl(kr)(XIM+Xl'_l)  +  ^yXhj(kr)(XIM-Xl'_l) 
2  2  k 
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Where  hj  (kr)  is  the  /th  order  spherical  Hankel  function  of  the  first  kind  and  the 
parameters  are,  respectively,  the  scattering  coefficients  for  the  TE  and  the  TM 

components. 

In  far-field  region,  where  kr  »  l ,  the  scattered  field  as  given  in  Eq.  (2)  can  be  simplified 
to: 
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In  Eq.  (4),  P'  is  the  associated  Legendre  function. 

From  Eq.  (3),  we  find  that  the  differential  scattering  cross  section,  which  is  defined  as  the 
optical  power  scattered  into  a  unit  solid  angle  and  normalized  by  the  intensity  of  the  input 
wave,  is  given  by: 


sc  _  dPsc 

dQ.  I  Inc  dkl 


•C£() 


■c£0  Eu 


cos-  </)+  5j  sin 


(5) 


The  total  scattering  cross  section  can  be  obtained  from  Eq.  (5)  by  integrating  the 
differential  cross-section  over  the  entire  4n  solid  angle  [16]: 


(Ttota/=-^E(2l  +  l)Rc[a/+J3/]  (6) 

k~  1=1 

From  Eq.  (6),  we  can  easily  find  the  total  scattering  cross-section  once  we  know  the 
coefficients  a/  and  fij  for  all  orders  of  partial  waves.  In  the  following  discussion,  we  outline 
a  transfer  matrix  method  developed  in  reference  [13]  and  apply  it  to  the  calculation  of  the 
scattering  coefficient  C//  and  /j/ .  To  begin  our  analysis,  it  suffices  to  consider  a  single  partial 
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wave  with  quantum  number  (l,  m),  since  the  incident  wave  and  the  scattered  wave,  as  given  in 
Eq.  (1)  and  (2),  are  formed  from  linearly  independent  combinations  of  components  with 
different  /  and  m  (m=  +  1 ).  For  the  / th  order  TE  or  TM  multipole  components  within  the  nth 
dielectric  layer,  we  can  express  the  corresponding  electromagnetic  field  as: 
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where  Zn  =  ^J/Jq  /  £§£n  is  the  material  impedance  and  kn  =  Je^co/  c  .  The  four  linear 
coefficients  An ,  Bn  C„,  and  Dn  are  constant  within  the  nth  layer.  Employing  the  continuity 
condition  of  Eq^E^^Hq  H ^  at  the  interface  between  two  adjacent  layers  and  the 
orthogonality  of  the  spherical  harmonics,  we  can  relate  the  linear  coefficients  An ,  Bn  Cn  and 
Dn  in  the  nth  layer  to  those  in  the  (n+  I  )th  layer  through  the  following  matrix  relations: 
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Applying  Eq  (8)  iteratively,  we  find  that  the  linear  coefficients  in  the  core  area  (Aco,  Bco, 
Cco,  Dco)  can  be  expressed  in  terms  of  those  in  the  outside  free  space  (Aol„,  Bou„  Cout,  Dolll)  as 


Bc, 

Ca 
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(9) 


In  Eq.  (9),  the  two  matrices  M™  and  MEE  are  obtained  from  multiplying  and  dividing  a 

series  of  two  by  two  matrices  given  in  Eq.  (8a)  and  (8b).  In  labeling  M  ™  and  MEE ,  we 
drop  the  quantum  number  in  since  the  matrices  are  identical  for  components  with  the  same  / 
but  different  m. 
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Combining  Eqs.  (1),  (2),  and  (9),  we  can  derive  a  simple  expression  for  scattering 
coefficients  a /  and  fij .  First,  we  notice  that  the  total  electromagnetic  field  outside  the  onion 

resonator  is  the  sum  of  the  incident  wave  and  the  scattered  wave  ( Etotai  =  Einc  +  Esc ),  which 
gives: 

°°  ;  J4;r(2/  + 1)  l~  or,  i  -  V  B,  \ 

Etot  =  Hi1  -  U,+^-hj)(XiM+XL_l)  +  -x(ji+^-hj)(XlM-XL_1)  (10) 

l=l  2  |_  2  k  2 

Meanwhile  according  to  Eq.  (7),  outside  of  the  onion  resonator,  the  electromagnetic  field 


for  the  component  with  quantum  number  (l,  m)  is  given  by: 

Elto?  ~  V  x  1 (Aoutj,  ( kr )  +  Bouthj  (kr))Xl  m  J  (TM)  (11a) 

~  (Coutj, (kr)  +  Douth\ (kr))Xl>m  (TE)  (lib) 

Comparing  Eqs.  ( 1  la)  and  ( 1  lb)  with  Eq.  ( 10),  we  find: 

A  out  ( Eout  =2 1  Pi ,  Cout  /  Dout  —  2/  a  i  (12) 

We  emphasize  that  this  relation  is  valid  for  all  multipole  orders.  Furthermore,  we  notice 
that  a  physical  solution  must  have  finite  value  at  the  origin,  which  leads  to 

BCo=(MfM)2'lAOut+(M?M)2aBOut=0  (13a) 

Dco  =  (MjE)2iCoul  +  (MjE)2aD0Ut  =  0  (13b) 

Combining  Eq.  ( 12)  and  Eq.  (13),  we  immediately  find 

Pi  /  2  =  —(M  ™  )  2  \  l(M  ™  )  2  2 

(14) 


/2  =  -(mJe)2  i  /(mJe)22 

Once  a /  and  /?/  are  obtained,  we  can  find  the  scattering  cross-section  of  the  /th  order 
partial  wave  and  the  total  scattering  cross-section  from  Eq.  (6). 

The  procedure  for  the  calculation  of  scattering  cross-section  as  outlined  above  applies  to 
cases  where  incident  wave  is  a  plane  wave.  In  experiments,  though,  the  incident  wave  is  likely 
to  have  a  Gaussian  beam  profile.  The  case  of  a  Gaussian  incident  wave  can  be  analyzed  using 
the  generalized  Mie  theory  [17,18],  where  the  effect  of  a  Gaussian  beam  can  be  included  by 

multiplying  the  scattering  coefficients  a )  m  and  Pj  m  by  a  beam  shape  factor  g’pjE  for  TE 

modes  or  g r{ljM  for  TM  modes.  In  this  paper,  we  limit  ourselves  to  the  simplest  case  where 
the  center  of  the  onion  resonator  coincides  with  the  center  of  the  waist  plane  of  a  fundamental 
Gaussian  beam.  In  this  case,  the  factors  and  g  are  reduced  to  the  same  factor  g / 
[18]: 


gl  =  exp- 


(/  +  1/2U]2 

2/rtT0 


(15) 


where  VT0  is  the  radius  of  the  minimal  spot  of  the  fundamental  Gaussian  beam,  A  is  the 
wavelength. 

To  describe  the  scattering  spectrum  of  the  onion  resonator,  it  is  often  convenient  to  use  a 
dimensionless  parameter,  the  total  scattering  efficiency,  which  is  defined  as  the  ratio  between 

the  total  scattering  cross  section  <7total  all(l  the  area  of  the  onion  resonator  cross-section: 
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In  this  definition,  R  is  the  radius  of  entire  onion  resonator.  Similarly,  we  can  define  the 
differential  scattering  efficiency  as  dQ^tt  / dQ.  =  d< 7SC  / jzR2  dQ. ,  with  dasc/dQ.  given  by 
Eq.  (5).  In  the  following  section,  all  numerical  results  are  given  in  terms  of  scattering 
efficiency  instead  of  scattering  cross-section. 

3.  Numerical  result  and  discussion 


In  this  section,  we  discuss  the  scattering  spectra  of  a  representative  Bragg  onion  resonator, 
which  is  composed  of  an  air  core  of  radius  rco  =  7 fan ,  followed  by  a  Si02  layer  and  Nr/a(i 
pairs  of  Si/Si02  Bragg  cladding  layers.  The  thickness  of  all  Si02  layers  is  Lsi(H  =  0.258/in , 
and  the  refractive  index  is  nsio^  =1.5  .  Correspondingly,  for  Si  layers,  we  have 

Lsi  =  0.1  ll/in  and  nsi  =  3.5  .  We  limit  ourselves  within  the  wavelength  range  of 
1.51-1.59  fan  .  The  bandgap  center  of  the  onion  cladding  is  located  at  1.55  flm  . 


Fig.  3.  Radial  dependence  of  the  electrical  field  of  TE  mode  and  the  magnetic  field  of  TM 
mode. 
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Fig.  4.  Spectrum  of  the  onion  resonator  modes  (calculated  with  Nciacj  =  4 ). 
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Before  calculating  the  scattering  spectrum  of  the  onion  resonators,  we  first  summarize 
some  important  properties  of  the  onion  resonator  modes.  As  mentioned  earlier,  the  onion 
resonator  modes  can  be  classified  into  core  modes  and  cladding  modes.  The  core  modes 
possess  relatively  small  angular  quantum  number  l  and  are  of  the  Fabry-Perot  type.  Their  field 
distribution  is  mostly  confined  within  the  onion  air  core.  In  contrast,  the  cladding  modes  have 
large  angular  quantum  number  l  and  are  similar  to  the  whispering-gallery  modes  in  disk 
resonators.  The  electromagnetic  fields  of  the  cladding  modes  are  concentrated  within  the 
cladding  layers.  In  Fig.  3,  we  show  the  radial  dependence  of  the  transverse  components 
(transverse  electrical  field  of  TE  mode  or  transverse  magnetic  field  of  TM  mode)  of  a 
cladding  mode  (TM42)  and  two  core  modes  (TEi0  and  TE24).  We  use  parameters  of  the  onion 
resonator  given  in  the  previous  paragraph  and  choose  N r/aij  =  4  .  The  features  in  Fig.  3 
clearly  demonstrate  the  different  localization  of  modal  energy  of  the  two  types  of  modes.  In 
Fig.  4,  we  plot  the  spectrum  of  the  core  modes  within  the  wavelength  range  of  1.51-1.59  fMn  . 
Similar  to  what  has  been  shown  in  Ref.  [14],  we  can  classify  the  core  modes  into  distinctive 
groups,  as  indicated  by  the  dashed  lines  in  Fig.  4. 


Fig.  5.  Spectrum  of  total  scattering  efficiency.  The  number  of  Bragg  pairs  is  Ncja(j  =  4  . 

The  scattering  spectrum  of  the  previously  discussed  onion  resonator  is  shown  in  Fig.  5.  In 
the  calculation,  we  assume  an  incident  plane  wave  and  take  N c/ad  =  4 .  In  the  figure,  the 
location  of  TE  and  TM  mode  frequencies  are  indicated  with  dashed  lines.  From  the  scattering 
spectrum,  it  is  clear  that  there  are  significant  differences  between  resonances  due  to  the 
cladding  modes  and  those  due  to  the  core  modes.  All  together.  Fig.  5  has  four  scattering  peaks 
that  correspond  respectively  to  the  TM42,  TM43,  TM44,  and  TM45  cladding  modes.  Their 
scattering  resonances  are  more  pronounced  in  amplitude  as  compared  with  those  of  the  core 
modes.  On  the  other  hand,  the  core  mode  resonances  exhibit  more  complicated  spectral 
features.  To  further  illustrate  the  difference  between  cladding  mode  resonances  and  core  mode 
resonances,  we  study  TEi0,  TE24,  TM42  and  TM43  modes  as  examples. 

We  use  the  scattering  formalism  developed  in  section  II  to  calculate  -Re(p;),  the  real  part 
of  the  /th  order  scattering  coefficient,  for  the  TM42  and  TM43  mode.  The  results  are  given  in 
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Fig.  6(a)  and  Fig.  6(b)  as  solid  dots.  Since  we  are  mainly  interested  in  the  resonant  part  of  the 
scattering  coefficient,  in  calculating  -Re(p,),  we  respectively  choose  /= 42  for  the  TM42  mode 
and  /= 43  for  the  TM43  mode.  As  in  the  case  of  whispering-gallery  modes  of  a  silica 
microsphere  [19],  the  resonant  scattering  coefficient  shown  in  Fig.  6(a)  and  Fig.  6(b)  can  be 
well  approximated  by  a  Lorentzian  line-shape  [19]: 

-  Re(«  or  /?)  = - - - -  (17) 

\  +  (Q-2AAl  Aq)2 


where  Q  and  Aq  are  respectively  the  quality  factor  and  the  wavelength  of  the  cavity  mode, 
and  AA,  is  defined  as  AA  =  A-A0.  In  Fig.  6(a)  and  Fig.  6(b),  we  plot  the  approximate 
Lorentzian  line-shape  as  given  by  Eq.  (17)  (shown  as  solid  lines).  We  emphasize  that  the 
Lorentzian  approximation  is  not  obtained  from  numerically  fitting  the  “exact'’  resonant 
scattering  coefficients.  Instead,  all  free  parameters  (Q  and  A0)  are  directly  calculated  from  the 
transfer  matrix  method  in  Ref  [13].  It  is  remarkable  to  notice  the  extremely  good  agreement 
between  the  two  results,  even  though  they  are  generated  from  two  different  approaches 
without  any  fitting. 


Fig.  6.  Comparison  of  scattering  resonances  numerically  obtained  using  the  Mie  scattering 
method  developed  in  section  2,  and  those  given  by  a  Lorentzian  approximation.  The  number  of 

Bragg  pair  is  Ncjaci  =  4  . 


The  total  scattering  efficiency  of  the  onion  resonator  has  contributions  from  both  a 
resonant  part  and  a  non-resonant  part.  As  a  first  order  approximation,  we  can  substitute  Eq. 


(17)  into  Eq.  (16),  add  a  constant  “background”  term  ,  and  write  the  total  scattering 
efficiency  near  resonance  as: 


f\scatt 

x^tot 


2/  +  1  2 
k2R2  l  +  (Q-2AA/Ao)2 


08) 


where  l  is  the  angular  quantum  number  of  the  resonant  onion  mode.  In  Eq.  ( 18),  only  Q^tt  is 

a  fitting  parameter,  whereas  Q  and  A0,  as  before,  are  all  directly  inferred  from  transfer  matrix 
calculations.  In  Fig.  6(c)  and  Fig.  6(d),  we  plot  the  total  scattering  efficiency  (calculated  from 
the  formalism  developed  in  section  II  and  represented  as  dots),  and  the  approximate  form 
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given  by  Eq.  ( 18)  (represented  as  solid  lines).  From  Fig.  6(c)  and  Fig.  6(d),  it  is  clear  that  the 
sharp  peaks  in  the  total  scattering  spectrum  can  be  entirely  attributed  to  the  presence  of  the 
cladding  modes,  and  we  can  accurately  extract  modal  frequencies  and  quality  factors  from  the 
scattering  spectrum. 

Significantly  different  from  those  of  the  cladding  modes,  the  spectral  characteristics  of  the 
core  modes  are  asymmetric  and  share  some  similarities  with  those  of  Fano  resonances  [20], 
From  the  scattering  resonances  shown  in  Fig.  5,  we  choose  to  study  two  representative  cases 
that  correspond  to  the  TEi0  and  TE24  mode.  In  Fig.  7(a)  and  Fig.  7(c),  we  give,  respectively, 
the  resonant  scattering  coefficients  -  Re[«;]  and  the  total  scattering  efficiencies  for  these  two 
modes.  In  the  figure,  the  corresponding  core  mode  wavelengths  are  indicated  using  dashed 
vertical  lines.  We  notice  that  for  the  core  modes,  the  wavelengths  with  maximal  or  minimal 
scattering  efficiencies  do  not  coincide  with  the  exact  modal  eigenwavelengths.  And  unlike 
cladding  mode  cases,  the  core  mode  scattering  spectra  do  not  have  the  usual  symmetrical 
Forentzian  line-shape.  Instead,  the  scattering  spectra  reach  the  minimal  (or  the  maximal) 
value  first  and  then  flip  quickly  to  the  opposite  limit.  We  also  notice  that  there  is  no  definite 
rule  regarding  whether  the  scattering  efficiency  first  reaches  a  dip  or  a  peak.  In  Fig.  7(b)  and 
Fig.  7(d),  we  plot  the  resonant  scattering  coefficients  and  the  total  scattering  efficiencies  for 
the  same  onion  resonator  but  with  Ndad=5.  It  is  interesting  to  notice  that  the  order  of  the 
appearance  of  the  peak  or  the  dip  changes  as  increases  from  4  to  5. 


4  layer 


5  layers 


4  layer 


5  layers 


Ainm) 


Fig.  7.  Scattering  resonances  of  the  TEio  and  TE24  modes  (with  N cjaij  =  4  and  5).  The 
dashed  lines  represent  the  resonant  scattering  coefficient  —  Re(dT;  )  (whose  values  are  shown 

on  the  left  y-axis),  and  the  solid  lines  represent  the  total  scattering  efficiency  (whose 

values  are  shown  on  the  right  y-axis).  The  dashed  vertical  lines  give  the  position  of  modal 
wavelength. 

In  Fig.  6,  we  have  demonstrated  that  for  the  cladding  modes,  we  can  extract  modal 
resonant  wavelengths  and  quality  factors  from  their  corresponding  scattering  resonances.  For 
the  core  modes,  even  though  their  scattering  spectra  is  no  longer  Forentzian,  we  can  still 
expect  a  definite  relationship  between  the  scattering  spectrum  and  the  corresponding  modal 
quality  factor  Q.  Quantitatively,  we  can  characterize  the  sharpness  of  the  core  mode  resonance 

using  the  modal  resonance  width  rlres  ,  which  is  defined  as  the  wavelength  difference  between 
the  peak  and  the  dip  in  the  resonant  scattering  coefficient  spectrum  (see  Fig.  7(a)).  Similarly, 
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we  can  define  the  modal  resonance  width  for  the  total  scattering  efficiencies.  In  Fig.  8, 
we  compare  A/Flres  ,  A  /  F^l  ,  and  quality  factors  of  TEi0  and  TE24  modes  as  a  function  of 
Nciad  ■  F°r  the  resonance  width  Tl°Js  calculated  from  the  total  scattering  efficiencies,  only 
results  with  Ncia(i  >  4  are  shown,  since  for  smaller  Ncia(/ ,  the  core  modes  have  low  quality 
factors  and  the  resonance  structures  due  to  different  modes  merge  together.  From  Fig.  8,  it  is 
clear  that  the  parameter  A/  Flrfs  provides  a  good  estimation  for  the  core  mode  quality  factors. 
As  a  result,  we  can  obtain  information  on  the  core  mode  wavelength  and  its  quality  factor 
from  the  scattering  spectrum  of  the  onion  resonator. 


Fig.  8.  Comparison  of  the  scattering  resonance  width  Fres  and  the  modal  quality  factor  as  a 

function  of  the  Bragg  pair  number.  F^es  and  F^s  represent,  respectively,  the  resonance 

width  of  the  scattering  coefficient  —  Re(flf)  and  the  total  scattering  efficiency  Qfgi'n  (see 
Fig.  7(a)). 


For  the  onion  resonator  shown  in  Fig.  1,  the  fabrication  process  can  introduce  a  finite 
amount  of  absorption  in  the  silicon  layers.  We  can  account  for  material  absorption  by 
introducing  an  imaginary  part  £t  to  the  Si  dielectric  constant.  Experimental  evidences  suggest 
that  Si  layer  absorption  is  of  the  order  of  dB/cm,  which  corresponds  to  a  value  of  £(~10"4. 
Since  the  cladding  modes  are  mostly  confined  in  the  cladding  region  and  the  core  modes 
mostly  within  the  air  core,  we  expect  that  the  presence  of  Si  absorption  will  have  bigger 
impact  on  the  cladding  modes.  In  Fig.  9,  we  compare  the  resonance  structure  of  TE24  and 
TM42  modes  for  three  different  values  of  £t :  0,  le-4  and  le-3.  As  expected,  the  influence  of 
absorption  loss  is  indeed  more  pronounced  for  the  cladding  mode  than  the  core  mode. 
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Fig.  9.  Effect  of  absorption  on  scattering  resonances.  The  number  of  Bragg  pairs  is 
N  clad  ~  4  ■ 

Finally  we  consider  two  issues  related  to  the  experimental  measurement  of  scattering 
spectrum.  As  discussed  previously,  the  input  wave  in  actual  measurements  is  likely  to  be  a 
Gaussian  beam,  which  can  be  accounted  through  multiplying  the  scattering  coefficients 
Of/  and  /?/  by  the  beam  shape  factor  g /  given  in  Eq.  (15).  In  Fig.  10(a),  we  compare  the 
scattering  spectrum  of  a  plane  wave  incident  on  an  onion  resonator  with  N c/ad  =  4  with  that 
of  a  Gaussian  incident  wave.  We  assume  that  the  center  of  the  minimal  beam  spot  of  the 
Gaussian  beam  coincides  with  the  onion  core  center,  and  use  a  minimal  beam  width  of  5  [dm 
in  the  calculation.  In  Fig.  10(a),  it  is  interesting  to  notice  that  the  Gaussian  incident  wave 
eliminates  the  scattering  resonances  due  to  the  cladding  modes  (i.e.,  TM42,  TM43,  TM44,  and 
TM45).  On  the  other  hand,  using  a  Gaussian  incident  beam  does  not  have  a  significant  impact 
on  the  core  mode  scattering  resonances.  Mathematically,  this  difference  between  the  core 
mode  and  the  cladding  mode  can  be  attributed  to  the  fact  that  the  beam  shape  factor  gh  as 
given  in  Eq.  (15),  clearly  favors  modes  with  smaller  angular  quantum  number  l  (i.e.,  core 
modes).  From  more  physical  point  of  view,  for  an  input  Gaussian  beam  with  a  smaller  beam 
width,  we  expect  a  lower  coupling  between  the  input  Gaussian  wave  and  the  onion  cladding 
modes,  since  the  cladding  modes  concentrate  in  the  region  relatively  further  from  the  center  of 
the  onion  resonator. 

In  experiments,  it  is  often  difficult  to  collect  light  scattered  by  the  onion  resonator  into  the 
entire  4n  solid  angle.  Instead,  we  consider  a  more  realistic  scenario  where  the  back-reflected 
light  within  the  range  of  0.9  n  <0<7t  is  collected.  The  corresponding  scattering  spectrum 
can  be  obtained  from  integrating  the  differential  scattering  cross-section  for  Q<(p<2n  and 
0.9/r  <0<7t .  In  Fig.  10(b),  we  show  two  back-scattering  spectra  that  correspond  to  the  two 
cases  investigated  in  Fig.  10(a)  (one  with  plane  incident  wave  and  the  other  one  with  a 
Gaussian  incident  wave).  As  in  Fig.  10(a),  we  find  that  using  the  Gaussian  incident  wave 
suppresses  the  resonances  due  to  the  cladding  modes.  Comparing  Fig.  10(a)  with  Fig.  10(b), 
we  notice  that  for  a  few  core  modes,  the  shapes  of  the  total  scattering  spectrum  and  the  back- 
scattering  spectrum  are  different.  However,  we  can  still  determine  the  wavelength  of  the  core 
modes  from  the  sharp  variations  of  the  back-scattering  spectrum. 
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total  scattering  efficiency 


back  scattering  efficiency  within  0:0.971-1 


X(jum) 

Fig.  lO.Comparison  between  plane  and  Gaussian  incident  wave.  In  both  (a)  and  (b),  the  upper 
line  is  under  the  assumption  the  incident  wave  is  a  plane  wave;  the  lower  one  is  under  the 
assumption  the  incident  wave  is  a  Gaussian  beam  with  the  waist  width  of  5  JLUn  with  the 
sphere  located  at  the  center  of  the  beam.  In  (a),  the  values  of  the  upper  and  lower  lines  are 
shown  on  the  left  and  right  y-axis  respectively.  The  number  of  Bragg  pair  is  Ncia(j  =  4  . 

4.  Conclusion 

We  combine  the  Mie  scattering  theory  and  a  transfer  matrix  method  to  analyze  light  scattering 
by  a  spherical  Bragg  onion  resonator.  We  present  and  analyze  numerical  results  for  several 
representative  cases.  We  find  the  scattering  resonances  due  to  the  core  modes  differ 
significantly  from  those  of  the  cladding  modes.  We  also  demonstrate  that  the  detailed  shape  of 
the  scattering  spectrum  depends  on  factors  such  as  the  incident  beam  profile  and  the  scheme 
for  the  collection  of  scattered  light.  Particularly  in  the  case  the  sphere  is  located  at  the  center 
of  the  beam,  the  cladding  modes  can  be  eliminated  by  shrinking  the  spot  size  of  the  Gaussian 
beam  while  the  influence  on  core  modes  is  not  significant.  We  conclude  that  the  wavelength 
and  the  quality  factor  of  the  onion  resonator  modes  can  be  obtained  from  the  position  and 
sharpness  of  the  resonances  in  the  scattering  spectrum. 
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